$dir = "/data2/bsi/RandD/bora/Exome_2samples/final_vcf";
$names=`ls $dir`;
#print $names."\n";
@names = split("\n",$names);
#print @names."\n";
$k = "/data2/bsi/RandD/bora/Exome_2samples/vcf-merge"; 
for($i=0;$i<@names;$i++)
{
	if($names[$i] !~ m/tbi/)
	{
		#$k = "tabix  -p vcf -f $dir/$names[$i]";
		#print $k."\n";
		#system($k);
		$k = $k." $dir/$names[$i]";
	}
}
$k = "$k | bgzip -c > /data2/bsi/RandD/bora/Exome_2samples/final_out.vcf.gz";
#print $k."\n";
system($k);
